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ABSTRACT 

In this paper, we study the heat transfer of MHD visco-elastic fluid flow through porous medium over an 
exponentially stretching sheet with viscous dissipation. By using suitable similarity transformations the governing partial 
differential equations are converted into non linear ordinary differential equations. These equations are analytically 
solved by using Homotopy Analysis Method (HAM). The convergence of the series solution which is obtained in the form 
of infinite series is discussed explicitly. The effect of various physical parameters like Magnetic parameter, Prandtl 
number, Eckert number, porosity parameter and visco-elastic parameter on velocity and temperature are discussed and 
displayed graphically. By increasing Prandtl number, the temperature decreases whereas it increases by increasing 
Eckert number. 
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1. INTRODUCTION 

The study of the boundary layer flow and heat transfer in viscoelastic fluids has ever increasing 
applications in engineering and manufacturing processes. Sakiadis [1] was the first to investigate boundary layer 
flow on continuous surfaces of finite length. Crane [2] extended this problem to a stretching sheet whose velocity is 
proportional to the distance from the slit which occurs in drawing of plastic sheet. P. S. Gupta and A. S. Gupta [3] 
analyzed momentum, heat and mass transfer of viscous fluid over a stretching sheet subject to suction/blowing. 
With this motivation Rajgopal et al. [4] examined the visco-elastic fluid flow over a stretching sheet. 

H. I. Anderson [5] studied the visco-elastic fluid flow over a stretching sheet in presence of transverse 
magnetic field. Later, M. I. Char [6] analyzed the heat transfer characteristics of an electrically conducting fluid in 
presence of transverse magnetic field. By taking into account the mass of the chemically reactive species, Sujith 
Kumar khan et al. [7] studied the heat and mass transfer of MHD visco-elastic fluid flow over a porous stretching 
sheet with viscous dissipation. P. S. Datti et al. [8] studied the effect of radiation of MHD visco-elastic fluid flow 
over a non iso- thermal stretching sheet with internal heat generation/absorption. Rafael Cortell [9] analyzed the 
flow and heat transfer of an electrically conducting visco-elastic fluid with viscous dissipation. M. S. Abel et al. 
[10] studied the visco-elastic fluid flow and heat transfer with viscous dissipation and non-uniform heat source. 
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All the above investigations are restricted to flow and heat transfer problems in non porous medium. The heat 
transfer in porous medium is carried out in many metallurgical processes. In these processes the filaments and continuous 
strips are cooled by drawing them through quiescent fluid. In the process of drawing, sometimes these strips are stretched. 
The rate of cooling can be controlled by drawing these strips in porous medium and the final product is obtained with 
desired accuracy. In this view, Vajravelu [11] studied the flow and heat transfer of viscous fluid in a saturated porous 
medium over an impermeable stretching sheet. Subhas et al. [12] studied the flow and heat transfer of visco-elastic fluid in 
porous medium over a stretching sheet. Mahantesh et al. [13] analyzed the heat transfer characteristics of visco-elastic fluid 
in porous medium over a stretching surface. 

Firstly, S. J. Liao [14] introduced Homotopy Analysis Method (HAM). Further it is developed and improved to 
solve non-linear problems. Subsequently, Liao [15] analyzed two branches of solutions to boundary layer flows over a 
stretched impermeable wall by using HAM method. Sajid et al. [16] applied Homotopy Analysis Method to discuss the 
analytic solution for MHD viscous flow due to a shrinking sheet. 

However, all these studies are confined to flow and heat transfer characteristics over a stretching sheet. 
By considering exponential stretching surface, Magyari et al. [17] described the heat and mass transfer characteristics of 
viscous fluid. Elbashbeshy [18] analyzed the heat transfer over exponential stretching surface with suction. Sajid et. al [19] 
studied flow and heat transfer of viscous fluid over an exponentially stretching sheet with radiation. Swati Mukhopadhyay 
[20] studied flow and heat transfer MHD viscous fluid over exponential stretching sheet in thermally stratified medium. 
Sujith Kumar Khan et al. [21] examined flow and heat transfer characteristics of visco-elastic fluid over an exponential 
stretching sheet. Hymavathi et al. [22] studied the momentum and heat transfer of MHD visco-elastic fluid over an 
exponential stretching sheet through porous medium by using Homotopy Analysis Method. 

The main objective of this paper is to study the flow and heat transfer characteristics of an electrically conducting 
visco-elastic fluid over an exponentially stretching sheet through porous medium with viscous dissipation by using 
Homotopy Analysis Method (HAM). The effect of various physical parameters like visco-elastic parameter (k[), porosity 
parameter (k 2 ), Magnetic parameter (M), Prandtl number (Pr) and Eckert number (Ec) on velocity and temperature are 
discussed through graphical representation. The results are compared with the available results in the literature and are seen 
in good agreement. 

2. MATHEMATICAL FORMULATION 


Consider an electrically conducting, steady, two dimensional boundary layer flow and heat transfer of viscoelastic 
fluid over an exponentially stretching sheet in porous medium. The flow is confined to y>0, x-axis is taken along the 
stretching surface in the direction of motion and y-axis is normal to it. By simultaneous application of equal and opposite 
forces along x-axis keeping the origin fixed the flow is generated as a consequence of exponential stretching of the sheet. 
A uniform magnetic field of strength B 0 is applied normal to the stretching surface. The flow and heat transfer 
characteristics under the boundary layer approximations are governed by the equations. 


du dv 
dx dy 


( 1 ) 
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Where u and v are the velocities in x and y directions respectively, T is the temperature of the fluid, p is the 
density of the fluid, k is the permeability coefficient of the porous medium, (7 is electrical conductivity of the fluid, C 
is the specific heat at constant pressure, 

The boundary conditions under the consideration are given by 

u = UJx) = U Q e x/l , v = 0, 7 = T w aty= 0 

u = 0, u y =0.T = T m + T Q e xl ' as y —> 

(4) 

Where U a ,T Q ,l are the reference velocity, the reference temperature and the reference length respectively. T rx is 
the temperature far away from the stretching sheet. 

Equation of continuity (1) is identically satisfied if we choose the stream function \j/ such that 

dif/ dy/ 

11=-^—, V =-— (5) 

dy c)x 


The momentum and energy equations can be transformed into the corresponding ordinary differential equations 
by introducing the similarity transformations: 


V = yy 




2 vl 


0 g x/2l 


(6) 


y/(x, y) = pvlU 0 f(x.T])e xl21 , 0(7]) = T (7) 

A W A oo 

Where TJ is the similarity variable / is the dimensionless stream function, 0(l]) be the dimensionless 
temperature 

The momentum and energy equations are transformed to 

r-2f' 2 +r-k^fr-\r v -lr 2 ]-(M + k 2 )f'=o (8) 

0 ''+ Pr(f0'-2f'0) + k 2 Pr Ecf' 2 +Pr Ecf" 2 = 0 (9) 

With the boundary conditions 
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f = 0, f' = 1, 9 = 1 at 7) = 0 

/' = 0, /" = 0, <9 = 0 astl^oo 


Where k, = 


a t U w 2vl 

- is the dimensionless visco-elastic parameter, k 2 = — - - 

pvl kU H 


is the porosity parameter. 


2 aB l l 

M = -— is the magnetic parameter, Ec 


pU 0 e 


x/l 


C P (T-TJ 


is the Eckert number, S = 


,a/2 UJ 


is the suction 


V 

parameter and rr = — is the Prandtl number. 

a 


3. HOMOTOPY ANALYSIS SOLUTION 


In this section, we employ HAM to solve the equation (8) and (9) subject to the boundary conditions (10). Let us 


choose the initial guesses f 0 ,6 0 as 

f 0 (tj) = l-e- rl (11) 

0 o (r/) = e’ 1 (12) 

The linear operator is selected as 

L(f) = f'"-f' (13) 

L(6) = 0" - 0 (14) 

Which have the following property 

L f [C x + C 2 e n + C 3 e~' 7 ] = 0 (15) 

L^CV+C^-^O (16) 


Where C ; (i = 1,2,3,4,5) are the arbitrary constants. 


If pe[ 0,11 is the embedding parameter, fl is the non zero auxiliary parameter and H (J/) is the auxiliary 
function, then we construct the following zero-th ordinary deformation equation: 

(1 - p)L(f (77 -p)-f 0 m = pfi,H X (p)Nffd 7 ; p)] (17) 

(\-p)L(6(T 1 -p)-e 0 (T 1 )) = pfi 2 H 1 (T 1 )N 2 [f(ri-p),e(T 1 -p)) (18) 

Subject to the boundary conditions 

f(0-p) = 0 f'(0\p) = 1 7 '( oo ; l ) = 0 
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We define non-linear operator as 
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For p= 0 and p = 1 , we have 
/(7/ ,0)= f 0 (T]), /(77,1) = /(7/) 

0(j],p)=0 Q ( 77 ) 6>(77,1) = < 9 ( 77 ) 


2ii 


(19) 


( 20 ) 


( 21 ) 


( 22 ) 


Thus as p increases from 0 to 1, f(t];p) varies from f (?] ;0) to f (t]) and 0(t];p) varies from 
0(f] ; 0) to 0(t]). Now, expanding f (t] , p), 0(t],p) in Taylor’s series with respect to p, we have 


/07 ./>) = /oW + Z/.W+' 

m=l (23) 

7(T7,/ ? ) = 7 0 (T7) + fX(7)7 m (24) 

m =1 


Where 
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■ V'f(thP) 
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1 d m m P ) 

m\ ()p m 


p =0 


(25) 


(26) 


If the initial approximations, auxiliary linear operators and non zero auxiliary parameters are chosen in such a way 
that the series (23) and (24) converges at p=l, then 


f(tj) = foW+Y.fjp) 

771=1 


(27) 


0(7) = 0 O (7) + IX (7) 

771=1 


(28) 
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Differentiating equations (17) and (18) ‘m’ times with respect to ‘p’, setting p=0 and ftnally diving with m!, we 
get the mth-order deformation equations are as follows: 

l, [/„ w - x. /„_, m=n ,h, wK w (29) 

4ie,";> - x,9, (7)1 = h 2 H 2 (ifiR‘.(.m oo) 

Subject to the boundary conditions 

/»(0) = 0 /I(0) = 0 /I(oo) = 0 


/>)=0 0 m (O) = O ^ m (oo) =0 


(31) 


m—1 m —1 m—1 i m—1 O m—1 

riw>=c, +Z/,/;.h -2i/,7;„-,-^[3Z/^.,/'--z/^-,/"--i/^-,/') 

i=0 i=0 j= 0 ^ z=0 ^ i=0 

-(fc 2 +M)/:_, 


(32) 
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(33) 


i=0 


i=0 
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We choose the auxiliary function as 
7 //^) = !, H 2 (T]) = \ 


(34) 


(35) 


If we let /*( 7]) , 0 m (T]) as the special functions of m lh order deformation equations, then the 


general solutions are given by 

/-(7) = /:(7) + C,+C 2 ^+C 3 e’ (36) 

0JrD = e; n (T]) + C A e^ (37) 

Where the integral constants Q(i =1,2,3,4,5) are determined using the boundary conditions (31). 

Now it is easy to solve the linear non-homogeneous equation (29) and (30) using MATHEMATICA software one 
after the other by considering m = 1, 2, 3... 

4. CONVERGENCE OF HAM SOLUTION 

By using Homotopy Analysis Method, the convergence region and rate of approximation of series solutions can 

, ti , and tl~. To find ti . andti 2 , tt curves are 

be adjusted and controlled by ustng the non-zero auxtliary parameters 1 2 1 2 

plotted in figure 1 and figure 2 for 20 ,h order of approximation. From figure 1 and figure 2 the range of admissible values 
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of tl | atld tl 0 are about — 0.91 < tl ^ <—0.32 and — 1.61</z, <—0.98. The convergence of HAM solution for 
different orders of approximations is given in Table.l. 




Figure 1: h-curve for f"(0) Figure 2: h-curve for 0 (0) 

5. RESULTS AND DISCUSSIONS 

In this paper, we analyze the heat transfer of MHD visco-elastic fluid flow in porous medium with viscous 
dissipation over an exponentially stretching sheet. The effects of visco-elastic parameter, porosity parameter, Magnetic 
parameter, Prandtl number and Eckert number on velocity and temperatures profiles are represented through graphs. 
Table.2 enlists the comparison of 070) for different values of visco-elastic parameter k, and are compared with the 
available results in the literature and are seen in good agreement. 

Table 1: Convergence of HAM Solution for Different 
Orders of Approximations when 

k } =0.1, k 2 =0.1, Ec = 0.1, Pr = 1.0 


Order 

-/"(0) 

-tf’(O) 

5 

1.424645 

1.203512 

10 

1.427125 

1.200031 

15 

1.427156 

1.199619 

20 

1.427156 

1.199521 

25 

1.427156 

1.199492 

30 

1.427156 

1.199483 

35 

1.427156 

1.199481 

40 

1.427156 

1.199481 


Table 2: Comparison of -0’(O) for Different Values 
of kj with k 2 = 0.5, Pr = 1.0, Ec = 0.25 


K| 

Mahentesh M (13) 

Present Study 

0.0 

1.12221 

1.122209 

0.1 

1.09631 

1.0963096 

0.2 

1.06551 

1.065513 
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n i 

Figure3: Velocity f' (q) for different Figure4: Temperature 0 (r|) for different 
Values of k, Values of k, 

The effect of visco-elastic parameter on velocity is shown in Figure 3. It can be shown that as viscoelastic 
parameter increases the velocity decreases i. e. the flow of the fluid retards with increasing of viscoelastic parameter. 
Figure 4 reveals the effect of visco-elastic parameter on temperature. It can be noticed that increase in viscoelastic 
parameter corresponds to an increase in temperature of the fluid. This is due to the fact that increase of visco-elastic normal 
stress causes the thickening of thermal boundary layer. 




Figure 5: Velocity f' (ij) for different Figure 6: Temperature 0 (ij) for different 
Values of k 2 Values of k 2 

Figure 5 is drawn for different values of porosity parameter k 2 on velocity. From figure 5 the velocity decreases 
with increasing values of porosity parameter, which leads to enhanced deceleration of the flow. Figure 6 indicates the 
effect of porosity parameter on temperature. It can be shown that by increasing the values of porosity parameter, the 
temperature of the fluid and thermal boundary layer thickness increases. 
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Figure 7: Velocity f' (ij) for different Figure 8: Temperature 0 (q) for different 
Values of M Values of M 

Figure 7 and 8 depicts the effect of magnetic parameter on velocity and temperature. It illustrates that by 
increasing the magnetic parameter the velocity decreases whereas temperature increases. With increase in values of M, rate 
of transport decreases due to the Lorentz force which opposes the motion of the fluid. 





Figure 9: Temperature 0 (q) for different 
Values of Pr 





Figure 10: Temperature 0 (q) for different 
Values of Ec 


The effect of Prandtl number and Eckert number is shown in Figure 9 & 10. By increasing Prandtl number the 
temperature decreases whereas it increases by increasing Eckert number. With increasing values of Prandtl number the 
thermal boundary layer thickness decreases so that the diffusion of heat is reduced. By increasing the values of Eckert 
number the thermal boundary layer thickness increases. This raise in temperature occurs as heat energy is stored in the 
fluid due to frictional heating. 

6 . CONCLUSIONS 


In this paper, we applied Homotopy Analysis Method to heat transfer of MHD visco-elastic fluid flow through 
porous medium over an exponentially stretching sheet with viscous dissipation. The effect of various physical parameters 
on velocity and temperature are discussed through graphs and the results are summarized as follows. 

• Visco-elastic parameter, porosity parameter and Magnetic parameter have same effect on temperature. i. e. 
temperature increases with increasing values of parameters. 
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• Diffusion of heat is reduced with increasing values of Prandtl number. 

• The thermal boundary layer thickness increases with increasing of Eckert number. 

• For effective cooling of stretching sheet visco-elastic liquids having low viscous dissipation must be chosen. 
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